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Abstract 


We use the Maple system to check the investigations of S. S. Gupta regarding the Smaran- 
dache consecutive and the reversed Smarandache sequences of triangular numbers [Smaran- 
dache Notions Journal, Vol. 14, 2004, pp. 366-368]. Furthermore, we extend previous 
investigations to the mirror and symmetric Smarandache sequences of triangular numbers. 
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The nth triangular number t,,, n € N, is defined by t, = )7}.,7 = n(n+1)/2. These numbers 
were first studied by the Pythagoreans. 


The first k terms of the triangular sequence {t,,}°°., are easily obtained in Maple: 


> t:=n->n*(nt+1)/2: 
> first := k -> seq(t(n) ,n=1..k): 


For example: 


> first (20); 
1,3, 6, 10, 15, 21, 28, 36, 45, 55, 66, 78, 91, 105, 120, 136, 153, 171, 190, 210 
In this short note we are interested in studying Smarandache sequences of triangular numbers 
with the help of the Maple system. 


To define the Smarandache sequences, it is convenient to introduce first the concatenation 
operation. Given two positive integer numbers n and m, the concatenation operation conc is 
defined in Maple by the following function: 


> conc := (n,m) -> n*10°length(m)+m: 


“Research Report CM04/D-05, June 2004. Submitted to the Smarandache Notions Journal. 
™MSc student of Informatics at the University of Craiova, Romania. Student at the University of Aveiro 
under the Socrates/Erasmus European programme, 2004. 


For example, 


> conc(12,345); 
12345 


Given a positive integer sequence {up }°°.,, we define the corresponding Smarandache Con- 
secutive Sequence {scs,}°°, recursively: 


SCS] = U1, 


SC8n = CONC(SCSn—1, Un) - 
In Maple we define: 


> scs_n := (u,n) -> if n=1 then u(1i) else conc(scs_n(u,n-1),u(n)) fi: 
> scs := (u,n) -> seq(scs_n(u,i),i=1..n): 


The standard Smarandache consecutive sequence, introduced by the the Romanian mathe- 
matician Florentin Smarandache, is obtained when one chooses un =n, Vn € N. The first 10 
terms are: 


> scs(n->n,10); 
1, 12, 123, 1234, 12345, 123456, 1234567, 12345678, 123456789, 12345678910 
Another example of a Smarandache consecutive sequence is the Smarandache consecutive 


sequence of triangular numbers. With our Maple definitions, the first 10 terms of such sequence 
are obtained with the following command: 


> scs(t,10); 


1, 13, 136, 13610, 1361015, 136101521, 13610152128, 
1361015212836, 136101521283645, 13610152128364555 


Sometimes, it is preferred to display Smarandache sequences in “triangular form”. 


> show := L -> map(i->print(i),L): 
> show([scs(t,10)]): 


1 
13 
136 
13610 
1361015 
136101521 
13610152128 
1361015212836 
136101521283645 
13610152128364555 


The Reversed Smarandache Sequence (rss) associated with a given sequence {u,}P2,, is 
defined recursively by 


TSS, = U1, 


PSS = CONC ig TAS Gh =7 ) 
In Maple we propose the following definitions: 


(u,n) -> if n=1 then u(1) else conc(u(n) ,rss_n(u,n-1)) fi: 


> rss_n := 
> rss := (u,n) -> seq(rss_n(u,i),i=1..n): 


The first terms of the reversed Smarandache sequence of triangular numbers are now easily 
obtained: 


> rss(t,10); 


1,31, 631, 10631, 1510631, 211510631, 28211510631, 
3628211510631, 453628211510631, 55453628211510631 


We define the Smarandache Mirror Sequence (sms) as follows: 


SMS] = U1, 


SMSn = CONC(CONC(Un, SMSn_—1), Un) 


> sms_n := (u,n) -> if n=1 then 

> u(1) 

> else 

> conc(conc(u(n) ,sms_n(u,n-1)) ,u(n)) 
> fi: 

> sms := (u,n) -> seq(sms_n(u,i),i=1..n): 


The first 10 terms of the Smarandache mirror sequence introduced by Smarandache are: 
> sms(n->n,10); 


1, 212, 32123, 4321234, 543212345, 65432123456, 7654321234567, 
876543212345678, 98765432123456789, 109876543212345678910 


We are interested in the Smarandache mirror sequence of triangular numbers. The first 10 
terms are: 


> sms(t,10); 


1, 313, 63136, 106313610, 1510631361015, 21151063136101521 
282115106313610152128, 3628211510631361015212836, 
45362821151063136101521283645, 554536282115106313610152128364555 


Finally, we define the Smarandache Symmetric Sequence (sss). For that we introduce the 
function “But Last Digit” (b1ld): 


> bld := n -> iquo(n,10): 
> bld(123) ; 
1 


If the integer number is a one-digit number, then function bld returns zero: 


> b1ld(3); 
0 


This is important: with our conc function, the concatenation of zero with a positive integer 
n gives n. 


> conc(bld(1) ,3); 
3 


The Smarandache Symmetric Sequence (sss) is now easily defined, appealing to the Smaran- 
dache consecutive, and reversed Smarandache sequences: 


$889n—1 = conc(bld(scsen—1),7S82n—1) ; 


S889n = CONC(SCS2n,1SS2n) , 


n €N. In Maple, we give the following definitions: 


> sss_n := (u,n) -> if type(n,odd) then 

> conc (bld(scs_n(u, (n+1)/2)),rss_n(u, (n+1)/2)) 
> else 

> conc(scs_n(u,n/2) ,rss_n(u,n/2) ) 

> fi: 

> sss := (u,n) -> seq(sss_n(u,i),i=1..n): 


The first terms of Smarandache’s symmetric sequence are 


> sss(n->n,10); 
1,11, 121, 1221, 12321, 123321, 1234321, 12344321, 123454321, 1234554321 


while the first 10 terms of the Smarandache symmetric sequence of triangular numbers are 


> sss(t,10); 
1,11, 131, 1331, 13631, 136631, 136110631, 1361010631, 1361011510631, 13610151510631 


One interesting question is to find prime numbers in the above defined Smarandache sequences 
of triangular numbers. We will restrict our search to the first 1000 terms of each sequence. 
All computations were done with Maple 9 running on a 2.00Ghz Pentium 4 with 256Mb RAM. 
We begin by collecting four lists with the first 1000 terms of the consecutive, reversed, mirror, 
and symmetric Smarandache sequences of triangular numbers: 


> st:=time(): Lscs1000:=[scs(t,1000)]: printf("%a seconds" ,round(time()-st)); 
20 seconds 


> st:=time(): Lrss1000:=[rss(t,1000)]: printf("%a seconds" ,round(time()-st)); 
75 seconds 


> st:=time(): Lsms1000:=[sms(t,1000)]: printf("%a seconds" ,round(time()-st)); 
212 seconds 


> st:=time(): Lsss1000:=[sss(t,1000)]: printf("%a seconds" ,round(time()-st)); 
26 seconds 


We note that scsj999 and rssjo999 are positive integer numbers with 5354 digits; 


> length(Lscs1000[1000]), length(Lrss1000[1000]) ; 
5354, 5354 


while smsjo99 and 88810999 have, respectively, 10707 and 4708 digits. 


> length(Lsms1000[1000]), length(Lsss1000[1000]) ; 
10707, 4708 


There exist two primes (13 and 136101521) among the first 1000 terms of the Smarandache 
consecutive sequence of triangular numbers; 


> st := time(): 

> select (isprime,Lscs1000) ; 

> printf("%a minutes",round((time()-st)/60)); 
(13, 136101521] 


9 minutes 


six primes among the first 1000 terms of the reversed Smarandache sequence of triangular 


numbers; 


> st := time(): 
> select (isprime,Lrss1000) ; 
> printf("%a minutes",round((time()-st)/60)) ; 
(31, 631, 10631, 55453628211510631, 786655453628211510631, 10591786655453628211510631] 


31 minutes 


only one prime (313) among the first 600 terms of the Smarandache mirror sequence of 
1 


triangular numbers; 
> length(Lsms1000[600]); # sms_{600} is a number with 5907 digits 
5907 


> st := time(): 

> select (isprime,Lsms1000[1..600]); 

> printf("%a minutes",round((time()-st)/60)); 
[313] 


‘Our computer runs low in memory when one tries the first 1000 terms of the Smarandache mirror sequence 
of triangular numbers. For this reason, we have considered here only the first 600 terms of the sequence. 


3 minutes 


and five primes among the first 1000 terms of the Smarandache symmetric sequence of trian- 
gular numbers (the fifth prime is an integer with 336 digits). 


> st := time(): 
> select (isprime,Lsss1000) ; 
> printf("4a minutes",round((time()-st)/60)); 


[11, 131, 136110631, 1361015212836455566789110512012010591786655453628211510631, 
1361015212836455566789110512013615317119021023125327630032535137840643546549652856159 
5630666703741780820861903946990103510811128117612251275132613781431148515401596165316 
5315961540148514311378132612751225117611281081103599094690386182078074170366663059556 
1528496465435406378351325300276253231210190171153136120105917866554536282115 10631 


19 minutes 


> length(%[5]); 
336 


How many primes are there in the above defined Smarandache sequences of triangular num- 
bers? This seems to be an open question. 

Another interesting question is to find triangular numbers in the Smarandache sequences of 
triangular numbers. We begin by defining in Maple the boolean function istriangular. 


> istriangular := n -> evalb(nops(select(i->evalb(whattype(i)=integer) , 
> [solve(t(k)=n)])) > 0): 


There exist two triangular numbers (1 and 136) among the first 1000 terms of the Smarandache 
consecutive sequence of triangular numbers; 


> st := time(): 
> select (istriangular,Lscs1000) ; 
> printf("%a seconds",round(time()-st)); 


[1, 136] 


6 seconds 


while the other Smarandache sequences of triangular numbers only show, among the first 
1000 terms, the trivial triangular number 1: 


> st := time(): 
> select (istriangular,Lrss1000) ; 
> printf("%a seconds",round(time()-st)); 


[1] 


6 seconds 


> st := time(): 
> select (istriangular,Lsms1000) ; 
> printf("%a seconds",round(time()-st)); 


[1] 
10 seconds 


> st := time(): 
> select (istriangular,Lsss1000) ; 
> printf("%a seconds" ,round(time()-st)); 


[1] 
6 seconds 


Does exist more triangular numbers in the Smarandache sequences of triangular numbers? 
This is, to the best of our knowledge, an open question needing further investigations. Since 
checking if a number is triangular is much faster than to check if a number is prime, we invite 
the readers to continue our search of triangular numbers for besides the 1000th term of the 
Smarandache sequences of triangular numbers. We look forward to readers discoveries. 
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